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■ This work continues our efforts to calibrate model surface brightness fluctuation 
| luminosities for the study of unresolved stellar populations, through the comparison 
. with data of Magellanic Cloud star clusters. We present here the relation between 

absolute K s -b&nd fluctuation magnitude and (V— /) integrated colour, using data from 
the 2MASS and DENIS surveys, and from the literature. We compare the s tar cluster 
sampl e with the sample of early- type galaxies and spiral bulges studied bv lLiu et alJ 
( 2002J . We find that intermediate-age to old star clusters lie along a linear correlation 
with the same slope, within the errors, of th at defined by the galaxies in the Mk b 
vs. (V — I) diagram. While the calibration bv lLiu et al.l was determined in the colour 
range 1.05 < (V - I c ) < 1.25, ou rs holds in the interval -5 > M y, > -9, 0.3 < (V - 1) 

K^J \ / W ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^M^^^^f"* ^^^^^^^ H ^^^^^^ H ^ H ^^^^^^^^^ H ^^^^^^ 

< 1.25 . This implies, according to lBruzual fc Ch ariot (200jj) and lMouhcine fc Lanconl 
( 2003) models, that the star clusters and the latest star formation bursts in the galaxies 
and bulges constitute an age sequence. At the same time, a slight offset between the 
^ ■ galaxies and the star clusters [the latter are ~ 0.7 mag fainter than the former at a 

given value of (V— /)], caused by the difference in metallicity of roughly a factor of two, 
confirms that the M k s vs. (V — I) plane may contribute to break the age-metallicity 
degeneracy in intermediate-age and old stellar populations. The confrontation between 
models and galaxy data also suggests that galaxies with K s fluctuation magnitudes 
that are brighter than predicted, given their (V — I) colour, might be explained in part 
by longer lifetimes of TP-AGB stars. A preliminary comparison between the H 2MASS 
data of the Magella nic star cluste r s and the sample of 47 early-type galaxies and spiral 
bulges observed by Uensen et alJ ((2003) through the F1QOW HST filter leads to the 
same basic conclusions: galaxies and star clusters lie along correlations with the same 
slope, and there is a slight offset between the star cluster sample and the galaxies, 
caused by their different metallicities. Magellanic star clusters are single populations, 
while galaxies are composite stellar systems; moreover, the objects analised live in 
different environments. Therefore, our findings mean that the relationship between 
fluctuation magnitudes in the near-IR, and (V — I) might be a fairly robust tool 
for the study of stellar population ages and metallicities, could provide additional 
constraints on star formation histories, and aid in the calibration of near-IR SBFs for 
cosmological distance measurements. 

Key words: astronomical data bases: miscellaneous — galaxies: star clusters - 
galaxies: stellar content — infrared: galaxies — infrared: stars — Magellanic Clouds 
- stars: AGB and post-AGB 
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1 INTRODUCTION 

Whereas surface brightness fluctuation 
ments are now well est ablished as ven 
indic a tors (e.g., |Tonrv 

120011: iMei et all 1200 it Uensen et all 120011 1200. 

potential for the study of unresolved stellar 



(SBF) measure- 
o werful distance 
Lhi^fe Graham! 



their 
popula- 



tions , although recognize d long since (|TQnrxA_ScI rngidaJ 
17983: iTonrv et~a7] S ILiu. Chariot, fc Ora.haml |200nt 
Blakeslee. Vazdekis. fc Aiharll200lt Ijensen et al.ll2003l) . has 

so far remained mostly unfulf ill ed. W e have recently 
jGonzalez. Liu, fc Bruzual 12004 120051: iMouhcine etaU 
l2005F ~started an effort to calibrate near-infrared (near- 
IR hereafter) SBFs, mostly for their use as probes of the 
characteristics of stellar populations. The idea is to com- 
pare SBFs derived from single stellar population (SSP) 
models with observed fluctuation luminosities of Magel- 
lanic Cloud (MC) star clusters. Our contribution has been 
to build "superclusters" with the data, in order to re- 
duce the stochastic effects produced by the inadequate 
representation, in single star clusters, of the luminos- 
ity fun ctions of stars evolving through short evolutionary 
phases llSantos fc FrogT^ll997l:lBruzuall2002l:ICervino et aU 

i 2002l: ICervino fc Valls-Gabaudl 120031: ICantiello et all 12005* 
Gonzalez. Liu, fc Bruzual 120041) ; i.e., the asymptotic giant 
branch (AGB) and upper red giant branch (RGB). Super- 
clusters, first introduced in lGonzalez. Liu, fc Bruz ual ( 2004, 
henceforth Paper I), are built by coadd ing MC clusters, in 
the compil ation of l van den BerghHl98ll). that have the same 
SWB class JSearle. Wilkinson, fc Bagnuololl980D . The SWB 
classification is based on two reddening-free parameters, de- 
rived from integrated ugvr photometr y of 61 rich star clus- 
ters in the Magellanic Clouds. Later, lElson fc Fail (Il985h 
assigned SWB classes to 147 more clusters using UBV pho- 
tometry. 

The SWB ranking constitutes a smooth, one- 
dimensional sequence that Searle et al. interpreted as 
one of increasing age and decreasing metallicity; this 
interpretation is supported by stellar photometry and 
spectroscopy in individual c l usters [e.g., references in 
Searle, Wilkin son, fc Bagnuolol dl98ot : iFrogel et al.1 lll990ft 1. 
The sequence discovered by Searle and collaborators was ar- 
bitrarily segmented by them to define a classification with 
seven types of clusters; however, facts like the apparently 
smooth progression in age and metallicity along the se- 
quence, and the similarity between the handful of type VII 
MC clusters and the old and metal-poor halo globular clus- 
ters of the Milky Way have translated into the customary 
adoption of the working hypothesis that clusters within each 
SWB type have approximately the same stellar population. 1 
In our case, as we have already explained in Paper I, binning 
the data into superclusters is the compromise we have found 
most convenient to adopt in order to try to circumvent the 
problem of small-number statistics posed by individual star 
clusters. We have devised eight superclusters, one for each of 



1 To give just a couple of examples. ICoher] Il982f) presents ages 
and metallicities of the seven SWB classes based on spectropho- 
tomctric d ata of a few stars w ithin only a couple of clusters of 
each class; IFrogel et alj ll99dl) do the same, based on near-IR 
photometry of 400 stars in and around 39 clusters. 



the seven different SWB classes, plus one "Pre-SWB-class" 
supercluster with the youngest objects in the sample. 

In the present paper, we focus on the relationship be- 
tween Mk s and the (V — I) colour. The merits of going 
to infrared wavelengths have been discussed at length [e.g., 
Liu et al. (2000); Paper I; Mouhcine et al. (2005)]. Sum- 
marising, the light of intermediate-age and old populations 
is dominated by bright and cool stars for which the spectral 
energy distributions peak in the near-IR, at the same time 
that dust extinction is reduced. On the other hand, there is 
a scarcity in the near-IR of both accurate empirical calibra- 
tions and self-consistent models, a fact that we hope we are 
contributing to remedy with our work. 

In the particular case o f the correlation between M k b 
and (V — I). \hiu et all (I2002T) confirmed a linear dependence 
of Mk s with (V — J) in a sample of 26 ellipticals, SOs, and 
spiral bulges in the Local Group, Fornax, Virgo, Eridanus, 
and Leo. The relation had already been hinted at by the data 
set of lJensen et al.l |l998), which covered only half the range 
in colour [~ 1.16 < (V — I) <~ 1.24]; it was more clearly 
shown bv lMei et alJ i200ll) . mainly t hrough the inclusion of 
data for NGC 4489 at (V - 1) < 1.05. lLiu et aU filled signifi- 
cantly the gap in colour with observations of Fornax objects. 
Probing the extent (in age, metallicity and, ultimately, en- 
vironment) to which this relation remains valid has obvious 
implications for an accurate calibration of near-IR SBFs as 
distance indicators, and hence for a p recise dete r minat ion of 
H through these measurements [see iLiu et al] (I2002T) for a 
detailed discussion]. 

There are also implications of this correlation from the 
point of view of stellar population studies. Since in general 
observed properties of stellar populations are luminosity- 
weighted, and given that after ~ 12 Myr the youngest 
stars are also the most luminous in the op tical and near-IR 
wavelengths, the tendency uncovered by ILiu et all J2002T) 
could be tracing the most recent burst of star formation 
in each of these systems. This is the conclusion reached 
in the cited work, after comparing their re sults with pre- 
dictions by ILiu. Chariot, fc Graham! i2000T) [based on the 



model s that would later be published by Bruzual & Chariot 
l200a)1. IWorthevHl993lJ). andlBlakeslee Vazdekis. fc Aihai 
i 200lli [based on the IVazdekis et al] (Il996ft modelsl. Al- 
though there are disagreements among such predictions (es- 
pecially in the derived metalliciti es obtained, respectively, 
from the lLiu. Chariot 1 fc Gra ham) models o n one hand, and 
from those of Blake slee. Vazdekis. fc AihaJ on the other), all 
of them point toward a large spread in age for the sample, 
from less than 5 Gyr to more than 12 Gyr. Furth e rmore , 
in all three sets of models confronted bv lLiu et al.l i2002l) , 
the age and metallicity sequences are not parallel in at least 
some regions of the Mk b vs. (V — 7) plane, opening the 
possibility of using it as a diagnostic for breaking the age- 
met allicity degeneracy. 

This study of the relationship between Mk b and (V — /) 
in the MC star clusters is comple mentary in two advanta- 
geous ways to the investigation bv lLiu et al.l i2002T) . Firstly, 
since the superclusters are approximately single age, sin- 
gle metallicity stellar populations, the star formation bursts 
that they represent are not masked by an underlying popu- 
lation, as it occurs in galaxies; hence, they likely constitute 
a better set for comparison with simple stellar population 
models. Secondly, the MC clusters span an even larger ex- 
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tent in age, i.e., from a few xlO 6 yr to ~ 10 10 yr, than that 
of the galaxies and bulges in the lLiu et alJ sample. Conse- 
quently, they cover a range roughly three times larger in 
M k s and four times more extended in (V — I). 

For this investigation, we have made use of 7f s -band 
data ret rieved from the Tw o Micron All Sky Survey 
C2MASS: ls"krutskie et aflll997b : J G „ (7 9 or I henceforth) 
data fro m the Deep Near-In frared Southern Sky Survey 
(DENIS: lEpchtein et al"lll997t l: and V data from different 
sources in the literature. This paper is organized as follows. 
In § [5] we give a brief summary of the data acquisition and 
characteristics, as well as of our own treatment of such data. 
In § |3| we present the ingredients of the stellar population 
synthesis models, and a comparison between the two sets of 
models used here. In §|3] we compare the theoretical predic- 
tions to the observations. Finally, in § [5J the results of the 
present work are discussed and summarised. 



2 OBSERVATIONAL DATA OF MAGELLANIC 
CLOUD STAR CLUSTERS 

2.1 The 2MASS data 

K s data for 1 91 M C clusters from the compilation of 
Ivan den Be rgh ( Il981_) and analised by Elson & Fall (1985, 
1988) were retrieved from the 2MASS archive. The data pro- 
cessing has been presented at length in Paper I. Succinctly, 
the fluctuation luminosity is the ratio of the second moment 
of the luminosity function (En^L 2 ) to its first moment (the 
integrated luminosity, T,rnLi), as expressed by the equation: 



Ew,L 2 



(i) 



Bright stars are the main contributors to the numerator, 
whereas faint stars contribute significantly to the denomi- 
nator. The second moment of the luminosity function was 
derived by measuring the flux of resolved, bright stars in the 
MC clusters, while the integrated luminosity was computed 
from the total light detected in the images, after removing 
the sky background emission. In order to reduce the stochas- 
tic errors produced by the small numbers of luminous and 
cool RGB and AGB stars in single star clusters, eight su- 
perclusters were assembled with the 2MASS data, one for 
each of the seven different S WB classes, plus one "pre-SWB- 
class" supercluster; this was accomplished by stacking in- 
dividual clusters with the same SWB type. 2 The mosaics 
were used to measure the integrated light of the superclus- 
ters. To derive the second moment of the luminosity func- 
tion, star lists for each one of the superclusters were inte- 
grated with entries from the 2MASS Point Source Catalog 
(PSC). The photometry of the point sources was performed 
by the 2MASS collaboration on individual cluster frames, 
previously to and independently from this work (i.e., we did 
not measure fluxes from point sources on the stacked super- 
cluster frames), following the standard procedure of profile- 
fitting plus curve-of-growth aperture correction. 3 Since all 

2 Individual cluster images were multiplicatively scaled to a com- 
mon photometric zero-point (determined by the 2MASS team) 
and dereddened; SMC clusters were geometrically magnified to 
place them at the distance of the LMC. 



near-IR zero-points were obtained in a uniform fashion, it 
is unlikely that they are an important source of systematic 
error. 

As explained in Paper I, the quoted fluctuation er- 
rors include stochastic variations produced by small-number 
statistics. These were calculated follow ing a statistical ap- 
proach introduced by iBuzzoni l)l989T) and ICervino et al. I 
(2002), based on the assumption that the variables involved 
have a Poissonian nature. In this framework, stochastic er- 

— 1/2 

rors scale as M tot , where M to t is the total mass of the 
stellar population. We refer the reader to the appendix for 
a short discussion on the subject of stochastic variations. 

Besides the problem of small-number statistics, assess- 
ing crowding is crucial for accurate SBF measurements of 
star clusters. Through the blending of sources, crowding 
can in principle make the numerator of equation Q larger 
and hence the SBF magnitude brighter. Another source of 
systematic error is the sky level, which impacts the de- 
nominator of equation Q Both crowding and sky determi- 
nation ha ve been lengthily addr essed in Paper I, and re- 
visited in lMouhcine et ail (|2005f ) to consider the corrected 
fluctuation measurements of su perclusters type I and II 
.Gonzalez. Liu, fc Bruzua3l2005h . In regard to crowding, we 
did not follow there the usual procedure of analising the 
effects of the addition of artificial stars on the measured 
fluctuations. Instead, we compared the fluctuations derived 
from annular regions of the superclusters (i.e., regions with 
diverse crowding properties), and checked whether the dif- 
ferent regions would be deemed cro wded by the criterion, 
developed by lAihar fc Tonrl (ll994Tl . that the two bright- 
est magnitudes of stars cover more than 2% of the area. 
We also tested our hypothesis that the PSC quality flags 
would be very helpful to eliminate blended sources by ob- 
taining, again in annuli, fluctuation values with and without 
stars with faulty photometry. While crowding affects prefer- 
entially the centers of clusters, errors in the sky subtraction 
will impact more the fainter, outer regions. Our analysis al- 
lowed us to determine that the circular regions within 1' of 
the center of the superclusters provide the better balance 
of uncertainties owing to crowding and sky subtraction, at 
the same time that they are less vulnerable to stochastic 
effects than smaller annular regions. Consequently, for the 
SBF measurements we use only integrated light and point 
sources within 1' (at the distance of the LMC) from the cen- 
tres of the superclusters. Just point sources with good pho- 
tometry, according to the PSC quality flags, 4 were included. 
We adopted distance moduli of (m — M) = 18.50 ± 0.13 
to the LMC and (m - M) Q = 18.99 ± 0.05 for the SMC 
.Ferrarese et al.ll2000D . 



2.2 The DENIS data 

The observations for DENIS were carried out between 1995 
and 2001, and i n particular the MC data were taken between 
1995 and 1998. ICioni et alJ <|2000T) give a good summary of 
the DENIS instrument, and data acquisition and character- 
istics. The instrument was mounted at the Cassegrain focus 
of the 1-m ESO telescope at La Silla, Chile, and could ob- 
tain simultaneously images at I g , J, and K s with 3 cam- 



' http:/ /www. ipac.caltech.edu/2mass/releases/allsky/doc/sec4_4. html 4 http:/ /www.ipac.caltech.edu/2mass/relcascs/allsky/doc/secl_6b.html#origphot 
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eras. These had, respectively, a Tektronix CCD with 1024 2 
pixels (each l"xl"), and two NICMOS infrared detectors 
with 256 2 pixels (each 3"x3", but the J and K s exposures 
were dithered to a l"pseudo-resolution). DENIS scanned the 
Southern sky in strips of 30° in declination and 12' in right 
ascension. Each strip comprises 180 12'xl2' images, with an 
overlap of 2' between every pair. The integration time at / is 
9 s; at J and K s , the integration time is 1 s, but each released 
image is made up of 9 individual microscanned exposures, 
for a total integration time also of 9 s. 

The nominal 5a limiting magnitude at I is 18 mag, and 
the typical size of a detected point source is smaller than 
2" FWHM. However, for this particular project we are con- 
cerned with the photometric accuracy at I achieved for ex- 
tended sources. From those clusters for which more than one 
/ calibrated exposure is available, we find the photometric 
error to be ~ 0.1 mag. 

Flattened and bias-subtracted / images that contain the 
Magellanic clusters in our sample were retrieved from the 
DENIS archive. The sample is the same one from which the 
Mk, measurements were obtained, with the exception of 
seven clusters, presented in Table 1; there are no calibrated 
DENIS data for six of them, while the other one (NGC 1777) 
has two bright foreground stars and we do not know whether 
the V measurement has been duly corrected. Instrumen- 
tal zero-points were obtained from the DENIS archive as 
well. Next, the /-band flux was derived for each cluster, in 
a diaphragm with the same size as the one used for the 
corresponding V-band me asurement. V magnitu des and di- 
aphragms were taken from lvan den Bergbl il98lfl , except for 
those clusters listed in Table 2. 5 It is worth noticing that, 
since we do not know the coordinates of the centres of the 
V observations, the derived (V — /) colour could be slightly 
biased to the red. 6 The sky emission for each cluster was de- 
termined from an annulus separated from the photometric 
diaphragm by a buffer area; the best sizes of both buffer and 
sky annuli were chosen after visually inspecting the images, 
with the purpose of excluding from the sky measurements 
bright foreground stars and, most importantly, residual clus- 
ter light. 

The individual cluster V and /-band flux values were 
then corrected for extinction as in Paper I, Table 1, and 
averaged to obtain the V and / fluxes for the superclus- 
ters, and their (V — I) colours. As before, when clusters do 
not have individually measured reddening, we have assumed 
E{B - V) = 0.075 for the LM C and E{B - V) = 0.037 for 
the SMC JSchleeel et alJll998l) . Once again, given that the 
extinction corrections of all the data (V, I, and K s ) were 
done consistently by us, this is probably not an important 
source of sys tematic error. Absolute zero poin ts were taken, 
for V from iBesselll <ll979T) . and for / from iFouaue et al] 
l|2000|b who determined it especially for DENIS. The av- 
erage (V — I) colours of the superclusters [and, for com- 
pleteness, also the Mk, values derived in Paper I and in 
iGonzalez. Liu, fc Bruzual i2005l) ] are presented in Table 3. 

5 In the case of the 2 clusters for which only F555W radial pro- 
files iMackev fc GilmordEoO^ were available, integrated magni- 
tudes were derived with diaphragms of 60". 

6 When we change the input centre for the / photometry -we do 
not have V images- between 5 and 10", the (V — I) colour of a 
single cluster shifts by about 0.1 mag. 



The quoted errors in (V — I) were derived as follows: the 
dispersion of / flux measurements was calculated for those 
clusters for which multiple images were available; for clus- 
ters with single images, the average dispersion was adopted; 
the error in / for the superclusters was found by adding the 
individual dispersions in quadrature, and then dividing by 
(TV — l) 1 / 2 , with TV the number of individual clusters in each 
supercluster; finally, the error in (V — /) was computed by 
assigning to the V value the same uncertainty as the one 
for /, and assuming a cor relation coefficient of 0.5 between 
V and / [see, for example. ICervino et al. I (I2002T) : Paper I]. 7 
We note, however, that while we find an average dispersion 
at / of ~ 0.1 mag for individual clusters, there might be 
systematics, like flatteni ng defects, tha t woul d increase this 
error somewhat. Indeed. IPaturel et al] <|2003T) give an aver- 
age uncertainty of ~ 0.2 mag at / for their measurements 
of galaxy magnitudes. We also remark that the dispersion 
of (V — I) colours among the individual constituents of each 
supercluster is typically ~ 0.2 - 0.3 mag. Although some 
fraction of this spread could be attributed to the fact that 
superclusters are not really single stellar populations, it is 
nevertheless consistent with the scatter among different low- 
mass realizations of a true SSP. F or example, using Monte 
Carlo simulations, iBruzuall <l2002i! finds 3-<r fluctuations in 
(V — K) of almost 2 mag for a lx lO 4 A/0 cluster at a given 
age. Moreover, the distribution of (V — I) colours of individ- 
ual clusters within each supercluster looks reasonably Gaus- 
sian for all except supercluster SWB VII. There are only 12 
star clusters type VII, and while the mode of the distribu- 
tion (3 clusters) lies around (V — /) = 1, there is a tail (2 
clusters) with (V — I) > 1.4. This tail could be produced 
by real population differences or, again, by small- number 
statistics. 



3 STELLAR POPULATION MODEL 
PREDICTIONS 

Central to this paper is the comparison between, on the one 
hand, the observed optical photometry and A" s -band fluctu- 
ation magnitudes of MC superclusters, and, on the other, the 
properties of synthetic single age, single metallicity stellar 
populations as predicted by, respectively, Bruzual & Chariot 
(2003, BC03) and Mouhcine & Langon (2003, ML03) stel- 
lar population synthesis models. BC03 models range from 
0.1 Myr to 17 Gyr in age, while those of ML03 go from 
12 Myr to 16 Gyr. Both sets of models span initial stellar 
metallicities of Z/Zq = 1/50 to Z/Zq = 2.5. Here, we de- 
scribe briefly the main ingredients of the stellar population 
synthesis models, referring to the quoted papers for more 
details. Note that other sets of theoretical predictions of 

7 The compilation bv lvan den B crgh ( 1981) docs not quote errors 
for the V photometry, but marks as uncertain (U — B) or (B — V) 
values for which any two observations dif fer by more than 0.1 mag. 
Amon g va n den Bergh's sources, only Ivan den Bergh &; Hager] 
il968t) and lAlcainol ll97Sfl list errors for V. Hence, the clusters 
in our work with published errors go from 18% for SWB class VI 
to 75% for SWB type VII. For supercluster type VII, we have 
calculated the uncertainty in (V — I) using the errors in these 
papers; it differs only by 0.01 mag from the one derived with the 
procedure described above. 
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near-IR _K" s -band surface brightness fluctuation magnitudes 
and optical colours are published; however, they generally 
use Bruzual & Chariot (2003) isochrones or one of their ear- 
lier versions (Liu et al. 2001; Mei et al. 2001), or they are 
computed for old stellar populations only, i.e., they do not 
cover the stellar population age range spanned by the MC 
superclusters (e.g., Worthey 1993a,b; Blakeslee et al. 2000; 
Cantiello et al. 2003). 

3.1 Single Stellar Population Models 

ML03 stellar population synthesis models were designed to 
reproduce the near-IR properties of both resolved and unre- 
solved stellar populations, with an emphasis on intermediate 
age stellar populations. The library of evolutionary tracks 
used by ML03 is based on the models of Bressan et al. (1993) 
and Fagotto et al. (1994 a,b,c). We will refer to these sets as 
the Padova tracks hereafter. The sets of tracks cover major 
stellar evolutionary phases: from the main sequence to the 
end of the early- AGB phase for low- and intermediate-mass 
stars, and to the central carbon ignition for massive stars. 
The Padova tracks do not extend to the end of the Ther- 
mally Pulsing Asymptotic Giant Branch phase (TP-AGB 
hereafter). The high luminosities and low effective tempera- 
tures of stars evolving through this phase make them among 
the main contributors to the integrated near-IR light of stel- 
lar systems within the age interval when these stars are 
alive (e.g., Frogel et al. 1990; Mouhcine & Langon 2002). 
The extension of these tracks to cover the TP-AGB phase 
is then needed. Until this happens, the evolution of low- 
and intermediate-mass stars through the TP-AGB evolu- 
tionary phase is followed using the so-called synthetic evo- 
lution modelling (e.g., Iben & Truran 1978; Renzini & Voli 
1981; see also Marigo et al. 2003 for another attempt to in- 
clude the TP-AGB in stellar population models). The com- 
plex interplay between different processes affecting stellar 
evolution through the TP-AGB phase are taken into account 
in the synthetic evolution models used by ML03. The prop- 
erties of TP-AGB stars are allowed to evolve according to 
semi-analytical prescriptions. These prescriptions take into 
account the effect of metallicity on the instantaneous prop- 
erties of TP-AGB stars. The evolution of TP-AGB stars is 
stopped at the end of the AGB phase, since their contribu- 
tion to the optical/near-IR light is almost negligible once 
they evolve beyond this phase. The models predict the ef- 
fective temperatures, bolometric luminosities, lifetimes, and 
relative lifetimes of the carbon-rich phase of TP-AGB stars. 

The BC03 models we will compare to the observed prop- 
erties of MC superclusters are those that these authors call 
the standard reference models. These models were built us- 
ing isochrones based on Padova stellar evolutionary tracks. 
To extend the Padova low- and intermediate-mass stars be- 
yond the early- AGB, BC03 adopt the effective temperature, 
bolometric luminosities, and lifetimes of TP-AGB stars from 
Vassiliadis & Wood (1993). The relative lifetimes in the 
carbon-rich phase for TP-AGB stars were taken from Groe- 
newegen & de Jong (1993), and Groenewegen et al. (1995). 

Stellar libraries are needed to assign spectral energy dis- 
tributions to the stars of a synthetic population. For stars 
evolving through phases other than the TP-AGB, ML03 
have used the theoretical stellar atmospheres of Kurucz (see 
Kurucz, 1979), Fluks et al. (1994) and Bessell et al. (1989, 



1991), as collected and re-calibrated by Lejeune et al. (1997, 
1998). BC03 use also the atmosphere compilations by Leje- 
une et al., as corrected by Westera et al. (2001, 2002). These 
stellar spectral libraries do not include spectra of TP-AGB 
stars. The spectrophotometric properties of TP-AGB stars 
are different from those of other cool and luminous stars. 
Their extended and cool atmospheres lead to the forma- 
tion of deep and specific spectral absorption bands. On the 
other hand, the photometric properties of oxygen-rich and 
carbon-rich TP-AGB stars are different, i.e., as an example, 
carbon stars show redder (H — K) colours than oxygen-rich 
TP-AGB stars at a given (J — H) colour. To account for TP- 
AGB star properties, ML03 have used the empirical library 
of average TP-AGB oxygen-rich and carbon-rich star spec- 
tra of Langon & Mouhcine (2002). The empirical average 
spectra of carbon stars from Langon & Mouhcine (2002) are 
used at all metallicities. It is extremely difficult, however, 
to estimate the effective temperature and the metallicity of 
TP-AGB oxygen-rich stars; instead, ML03 use a metallicity- 
dependent effective temperature scale (Bessell et al. 1991). 
To account for the spectral properties of carbon-rich TP- 
AGB stars, BC03 have constructed period-averaged spectra 
using solar metallicity model atmospheres for carbon stars 
from Hofner et al. (2000). These spectra were used to rep- 
resent carbon stars of all metallicities in BC03 models. 

i^s-band fluctuation magnitudes were computed from 
the models , as has been amply des cribed in Paper I and sum- 
marised in lMouhcine et alJ (§2005) . In this case, the second 
moment of the luminosity function is calculated by summing 
L 2 Ks over all stellar types in the models; the integrated lumi- 
nosity is obtained by summing the Lk b of all types of star. 
Correspondingly, model (V — /) colours were constructed by 
summing the V and I fluxes of all the stars in the isochrones. 

3.2 Comparison between models 

Fig.Qpresents a comparison between the evolution of single- 
age, single-metallicity stellar population properties in the 
M Ks vs. (V - I) diagram as predicted by BC03 and ML03 
stellar population synthesis models, respectively, for all the 
metallicities considered. The figure shows that both sets of 
models predict a qualitatively similar evolution over the age 
range comprised by them; i.e., the Jf s -band surface bright- 
ness fluctuation magnitudes get fainter as the (V — I) colour 
gets redder. For both sets of stellar population model predic- 
tions, the (V — I) colour increases gradually to redder values 
as the stellar populations age. On the other hand, it is pre- 
dicted by both BC03 and ML03 models that at a fixed age, 
single age, single metallicity stellar populations show red- 
der (V — I) colours at higher metallicity. This is because, at 
fixed initial stellar mass, lowering metallicity causes stars to 
evolve at higher effective temperatures. In view of the mono- 
tonic and smooth evolution of the (V — I) colour, and given 
its weaker sensitivity (compared to the near-IR wavelength 
range) to the presence of very cool and luminous stars, this 
colour index can be regarded as a primary age indicator at 
fixed metallicity. 

For the stellar populations dominated by red supergiant 
stars, i.e., younger than a few x lOMyr, the _ff s -band SBF 
magnitude gets drastically fainter with age; on the other 
hand, the (V — I) changes by a modest factor around ~ 0.3 
mag. This trend is produced by the combination of two facts: 
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(1) the masses, and hence the luminosities, of the red super- 
giant stars that drive the fluctuation signal change signifi- 
cantly, while (2) the combined mass of main sequence stars, 
that determine the integrated optical properties, stays al- 
most constant. 

When red supergiant stars disappear from a stellar pop- 
ulation, AGB stars drive the evolution of the near-IR prop- 
erties up to 1.5-2 Gyr. For populations dominated by short- 
lived massive AGB stars, i.e., younger than ~ 200 Myr, 
the predicted evolution in the M k b vs. (V — /) diagram 
is complex. When the first (massive) AGB stars emerge in 
the stellar population, a brightening of the AVband sur- 
face brightness fluctuation magnitude is predicted, at almost 
fixed (V — I) colour. This is because of the overluminosity 
produced by the envelope burning that affects AGB stars 
with large initial stellar masses, i.e., Minit > 3.5 — 4Mq (see 
e.g., Mouhcine & Lancon 2002 for more details on the effects 
of envelope burning on intermediate-age stellar population 
properties). The observed counterparts of single age, sin- 
gle metallicity stellar populations within this age range are 
expected to cluster at the same location in the Mk s vs. 
(V — I) diagram, i.e., in the region around (V — /) ~ 0.4 
and M k s ~ —7.5. During this age interval, the models have 
no ability for stellar population age-dating, and/or metallic- 
ity estimate. 

For stellar populations older than ~ 300 Myr, in which 
AGB stars are not affected by envelope burning, both sets of 
models predict a monotonic dimming of the A" s -band surface 
brightness fluctuation magnitudes as the (V — I) colour gets 
redder. This evolutionary pattern continues for older stel- 
lar populations, i.e., older than 2-3 Gyr, when RGB stars 
drive their near-IR properties. This is because of the evolu- 
tion of late-type giant star content. As a stellar population 
ages, the mass of the stars fueling the evolution of near-IR 
properties, i.e., AGB stars for ages younger than ~ 1.5 Gyr, 
and red giant stars for older ages, decreases. Consequently, 
at a fixed metallicity, the average luminosity of these stars 
decreases with stellar population age. Conversely, at a fixed 
age, M k s magnitudes get brighter and the (V — I) colour 
gets redder as the stellar population metallicity increases. 
Thus, both the BC03 and ML03 sets of models predict that, 
as the metallicity increases, populations with the same age 
move to the upper right in the Mk s vs. (V — I) diagram. 

Despite the qualitative agreement between the BC03 
and ML03 sets of theoretical predictions, differences between 
the two are apparent. At a given stellar metallicity, the evo- 
lutionary track in the Mk s vs. (V — I) diagram predicted by 
ML03 models is systematically redder than the one predicted 
by BC03 models. For stellar populations with Z= 0.0004, 
0.008, and 0.02, and between ages 300 Myr and 1.5 Gyr, 
the models based on ML03 isochrones predict, at similar 
(V — I) colours, bright er iG-band SBF magnitudes [see also 
iMouhcine et alJ i2005Tt l . Observationally, this means that for 
a given (V — I) colour, models based on BC03 isochrones 
will attribute a higher metallicity to a star cluster with a 
certain M k s magnitude. The A" 3 -band SBF magnitude in 
the models based on ML03 isochrones is more sensitive to 
the presence of AGB stars. This is due to the fact that the 
AGB lifetimes used in the ML03 stellar population synthesis 
models are larger, thus increasing the contribution of these 
stars to the Kg-band light budget. 



4 MODELS VERSUS OBSERVATIONAL DATA 

Figures|2jand|31show the comparison of Mk 3 versus (V — I g ) 
colour of MC superclusters with both sets of models. BC03 
models are presented in the left panels, while ML03 ones 
are displayed on the right. Fig. [5] presents models with Z = 
0.0004, Z = 0.004, and Z = 0.008. Even though the super- 
clusters all have Z < 0.01, for completeness Fig.[3]illustrates 
the evolution of models with solar metallicity and Z — 0.05; 
the models with Z — 0.008 are also shown in Fig. |3] for 
comparison purposes. The solid dots represent the MC su- 
per cluster data, and the rectangle marks the general locus of 
the lLiu et alJ (I2002T) galaxy sample. T he observed (V — I c ) 
galaxy colours reported by iLiu et al] were transformed to 
(V — I g ), assuming that the mean stellar metallicities of the 
sample galaxies cover the range 0.008 < Z < 0.05, by the 
following equation: 

(V - I g ) = (1.050 ± 0.006)(V - Ic) - (0.005 ± 0.002) (2) 

Focusing first on the data and on Fig. we notice that 
the youngest Pre-SWB and SWB I superclusters, respec- 
tively at (V - I) = 0.34 and (V - I) = 0.61, are observed 
with significantly redder optical colours than predicted by 
the models with Z — 0.008, the o nes closest to their metal- 
licity of Z = 0.01 JCohenl 11982^ Given their young ages, 
though, this is not surprising; Tcharlot fc Fall ( 2000) offer 
the prescription that populations younger than ~ 10 7 years 
suffer from about three times more reddenin g than later in 
their lifetimes. 8 Concordantly, for example, iGrebel fc Chul 
(2000) have measured a total colour excess (including fore- 
ground extinction) E(B -V) = 0.28±0.05 for Hodge 301, a 
relatively old cluster in 30 Dor for which these authors also 
derive an age of 20-25 Myr. This is the mean age of the star 
clusters that compose the supercl uster type I (Paper I), if 
one adopts the age calibration by lElson fc Fall (1985J). In- 
terestingly, if one dereddens the SWB I supercluster by an 
extra 0.20 mag i n E(B — V), o r the d ifference between the 
measurement bv lGrebel fc Chul feOOOD and the average 0.08 
mag that was adopted for the LMC dSchlegel et alJ ll998). 
the data point falls exactly on the models with Z = 0.008. 
We f ollow a s i milar proc edure for the Pre-SW B supercluster. 
Both Parker (1993) and Dick el et al] ljl994fl have measured 
an average reddening of E(B — V) ~ 0.43 towards the whole 
of the 30 Dor region, while ISelman et al] lll999fl have de- 
termined a radially dependent reddening in the direction of 
R136 that reaches E(B — V) > 0.5 at the center of the clus- 
ter and declines to E(B — V) ~ 0.3 at r = 1'. Dereddening 
the youngest, Pre-SWB supercluster by E(B — V) = 0.35 
above the average of the LMC places it right where the BC03 
models predict a population younger than 5 Myr should be 
in Figures [5] and [3] It is possible that the fluctuation magni- 
tudes derived for the Pre-SWB and t ype I superclusters su f- 
fer from a certain degree of crowding l|Mouhcine et al .1200511 . 
If this is the case, the correct values would be fainter than 
the reported ones, but would still lie along the same model 
sequence, at ages, respectively, slightly younger and slightly 
older. For the rest of the MC supercluster data points we 
attempt no further correction. 

8 Exactly the same relation between age and r edden ing was found 
observationally by van den Bergh & Hagcn (1968) for the MC 
clusters. 
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Figure 1. Comparison of model M Ks versus (V ~ I g ) colour for Z = 0.0004 (top left), Z = 0.004 (top middle), Z = 0.008 (top right), 
Z = 0.02 (bottom left), and Z = 0.05 (bottom middle). Solid lines: Bruzual & Chariot (2003, BC03); dotted lines: Mouhcine & Lancon 
(2003, ML03). For Z = 0.0004 and Z = 0.004, evolution of models is shown only between ~ 300 Myr and ~ 16 Gyr. In the remaining 
panels, BC03 and ML03 models start, respectively, at 2.4 Myr and 12 Myr. From left to right, tickmarks indicate 5 and 10 Myr (only 
for BC03, Z > 0.008 models), 100 and 200 Myr (for BC03 and ML03 models with Z > 0.008), 0.4, 1, 4, and 12 Gyr (all models). 



Table 1. Clusters in Gonzalez et al. (2004) not included in this paper 



Supercluster 


Name 


Reason 


I 


L 51 


No available calibrated DENIS data 




SL 477 


No available calibrated DENIS data 




NGC 1951 


No available calibrated DENIS data 




NGC 1986 


No available calibrated DENIS data 


Ill 


NGC 1953 


No available calibrated DENIS data 


V 


NGC 1777 


Two bright foreground stars 




NGC 2193 


No available calibrated DENIS data 
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Table 2. V photometry sources different from Ivan den Berghl Jl98 ill 



Supercluster Name Source 



Pre-SWB . . NGC 1727 Bica et al. f!996l 

NGC 1936 (IC 2127) Bica et al. (19961 



NGC 2014 Bica et al. ("19961 

NGC 2074 Bica et al. (19961 

NGC 2001 Bica et al. (19961 

I NGC 299 Alcaino (19781 

IV SL 663 Mackev fc Gilmore (20031 ° 

V SL 363 Bernard k. Bigav (19741 : Bernard (19751 ; Bica et al. (19961 

SL 556 Bica et al. (19961 

SL 855 Mackev fc Gilmore (20031 ° 

VI SL 842 Bica et al. (19961 

ESO 121-SC03 Mateo. Hodge, Schommer (19861 

VII NGC 1786 van den Bergh (19811 and Bica et al. (19961 b 



a F555W magnitudes we re first obtained by integrating over the radial luminosi t y pro files published by 
iMackev &; Gilmord 120031) . then transformed to V with the rela tions of iDolphinl l2000l) (wi thout charge 
trans fer terms), assuming (B-V) colours from lElson &: Falll ( 1985) appropriate for the s types felson fc Falll 
ll98Sft of the clusters. 

b Ivan den Be Tghl <198lD lists V = 10.88 and states that only observations corrected f or the contribution of 
a foreground star have been considered, but neglects giving a diaphragm. iBica et alj Jl996l) register both a 
V mag of 10.88 and a diaphragm of 60". 



Table 3. M Ks and (V - I) values 



Supercluster Mk s (V — I) 



pre -7.70±0.40 0.34 ± 0.14 

I -8.85±0.12 0.61 ± 0.08 

II -7.84±0.28 0.48 ± 0.10 

III -7.45±0.24 0.47 ± 0.06 

IV -7.51±0.18 0.54 ± 0.05 

V -6.69±0.20 0.78 ±0.11 

VI -6.21±0.24 1.02 ± 0.07 

VII -4.92±0.38 1.06 ± 0.13 



The evolution predicted by the models of the K s -baiid 
SBF magnitudes as a function of the (V — I) colour agrees 
remarkably well with the observed sequence defined by the 
MC superclusters. Superclusters of SWB types II, III, and 
IV have ages, respectively, of 35, 105, and 320 Myr, us ing the 
SWB class-age transformation of iFroeel et alJ il990fl . Their 
observed loci concur nicely with the predicted location of 
stellar populations with ages between a few x ~ 10 Myr 
and ~ 300 Myr, in the region at (V — I) ~ 0.4 mag and 
M k b ~ —7.5 mag. SWB types V, VI, and VII correspond 
to ages of ~ 1 Gyr, 3 Gyr, and 9 Gyr, respectively. For older 
stellar populations, the predicted monotonic fading of the 
JG-band SBF magnitudes as the (V — I) colour gets redder 
agrees with the observed properties of intermediate-age and 
old MC superclusters. 

However, one cannot fail to notice that the MC super- 
clusters older than ~ 3 Gyr (i.e, with SWB types VI and 
VII) are not very well matched by the models. Given their 
respective ages and metallicities, 9 the type VI superclus- 

9 For these SWB types, 0.002 < Z < 0.0008 <Frogel et alll99Cft 



ter mostly seems brighter in Mk s , while type VII mainly 
appears redder in (V — I) than the model predictions. We 
remind the reader here that the measured (V — I) colours 
could be slightly biased to the red as a consequence of our 
ignorance of the exact telescope pointings used for the ac- 
quisition of the V data (see § 12.21 . We have discussed in 
the same section that the colour distribution of the individ- 
ual SWB VII clusters shows evidence of slight population 
differences, although we could also be seeing the effects of 
small- number statistics (12 clusters). 

As advanced in § Q it is readily apparent in Figures 
|5| and |3| that the MC superclusters have greatly increased 
the range in Mk b versus (V — I) that is now accessible for 
exploration. We have plotted in both pan els of F ig. [3 as a 
solid straight line, the linear fit found bv lLiu et alJ (2002) 
to their early-type galaxy and spiral bulge data set. After 
transforming (V — I c ) to (V — I g ) with eq. |2] this relation 
becomes 

M K , = (-5.84 ± 0.04) + (3.4 ± 0.8)[(V - I g ) - 1-20]. (3) 

It is remarkable that, with the exception of the ex- 
tremely young supercluster Pre-SWB, and the old, metal- 
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Figure 2. Comparison of Mk s versus (V — I g ) colour of Magellanic Cloud superclusters with low-Z models by BC03, left panel), and by 
ML03 right panel). Dotted line: model with Z = 0.0004; dashed line: Z = 0.004; solid line: Z = 0.008. Tickmarks as in Fig.0 The thick 
arrows deredden the Pre-SWB and SWB I superclusters by E(B — V) = 0.2 above the average of the L MC. The thin arr ow dereddens the 
Pre-SWB supercluster by an additional E(B — V) = 0.15 (see text). The rectangle marks locus of the lLiu et alJ 12002?) galaxy sample. 
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Figure 3. Comparison of M k 3 versus (V — I g ) colour of Magellanic Cloud superclusters with high-Z models by BC03 (left panel), and 
by ML03 (right panel). Solid line: model with Z = 0.008; dotted line: Z = 0.02; dashed line: Z = 0.05. Evolution before ~ 300 Myr is 
now shown for all models. Symbols as in Figures and [5] except that BC03 model with Z = 0.02 is graphed without tickmarks, since 
it runs so closely to the Z = 0.008 model. The straight solid line is best fit to lLiu et al.l 120021) sample of early-type galaxies and spiral 
bulges. The straight dotted-dashed line is our fit to superclusters I through VI (see text). 



poor supercluster SWB class VII, the remaining MC super- 
cluste rs lie roughly along the correlation found bv lLiu et alJ 
(2002). In fact, a fit to superclusters I (corrected for extra 
extinction as explained in Q through VI yields: 



This is also shown in Fig. [3] as a straight dotted-dashed 
line. Within the errors, the relations for, respectively, MC su- 
perclusters and galaxies have the same slope, 10 at the same 



M Ks = (-5.1 ± 0.4) + (4.0 ± 0.6)[(V - /, 



1.20]. 



(4) 



10 If in fact the fluctuation magnitude of supercluster type I is 
affected by crowding and its true value is a few tenths of a mag 
dimmer, the match between the two slopes would be even better. 
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time that the MC superclusters have a systematically lower 
M k b magnitude, at a fixed (V — I) colour, than what is 
predicted for early-type galaxies. The observed offset can be 
explained, according to both the BC03 and ML03 sets of 
models, as a metallicity effect. The metallicity of the MC 
superclusters is roughly (bar only S WB class VI) half that 
of the sample of early-type galaxies and spiral bulges. For a 
given (V — I) colour, the metal-poor simple stellar popula- 
tions are on average older than the metal-rich ones, having 
then fainter M k b as a consequence of both lower metallici- 
ties and older ages. 

We refer the reader to Paper I for a detailed discus- 
sion of the reasons for the monotonic decline of the K s - 
band SBF brightness with age displayed by intermediate- 
age and old populations. When compared to both sets of 
models, the result obtained here fo r the MC s uperc lusters 
backs the interpretation offered in iLiu et al.1 (l2002l) . that 
their galaxy sample constitutes a sequence. Namely, that 
objects with brighter SBFs have a more recent latest star 
formation burst (or a more extended star formation his- 
tory), and that the most recent bursts in every object have 
all occurred at roughly constant metallicity. This is so be- 
cause their sample lies on a straight line, parallel to the 
age sequence of their single metallicity models. Simultane- 
ously, the fact that there exists a discernable offset between 
the galaxies and the MC superclusters confirms that the 
M k b vs. (V — I) plane may contribute to the decoupling 
of age and metallicity effects in intermediate-age and old 
stellar systems. Particularly worthy of notice is the observa- 
tion that the two samples offer consistent results, considering 
that galaxies are composite stellar systems, while the super- 
clusters are approximately single stellar populations. The 
most likely explanation is that, being mostly probes of the 
brightest stars in a population at a given wavelength, SBFs 
are relatively insensitive to an underlying older population 
in composite systems. 

It is also interesting to remark here that, while the 
BC03 models attribute a higher metallicity to NGC 1419 
and NGC 1389, the two Forn ax objects left o ut of their 
fit (through sigma-clipping) bv lLiu et all J2002I) . the ML03 
model with solar metallicity (Fig. [3] right panel, dotted line) 
shows an excursion to brighter SBF magnitudes at exactly 
the right place to ascribe to them the same metallicity as 
to other galaxies that do lie on the linear fit. This is a con- 
sequence of the longer lifetimes of TP-AGB stars in ML03 
models. Hence, longer TP-AGB lifetimes could partly ex- 
plain the observations of galaxies with M k b brighter than 
predicted, given their (V — /) colour. 

Last but not least, we want to mention here that the 
largest sample to date of gala xies put together fo r SBF 
studies has been presented by Ijensen et al.l ((20033 • The 
65 objects in this sample, however, were observed with 
the Near-IR Camera and Multi-Object Spectrometer (NIC- 
MOS) camera 2 (NIC2) on board the Hubble Space Tele- 
scope (HST), using the F160W filter. The photometric 
transformations between the NIC2 and the commonly used 
near-IR ground-based filters is particularly difficult, for sev- 
eral reasons: the very deep molecular absorption bands in 
the star themselves; the significant differences between the 
HST and the ground-based filters; and the fact that tel- 
luric absorption features are very deep in ground-based 
observations but completely absent in data obtained from 



space dStephens et al-lboOOT) . Dealing with fluctuation mag- 
nitudes has an added degree of complexity, and conflict- 
ing statements can be found in the literature with regard 
to what colour to use to determine the correct transfor- 
mation coefficients [i.e . , the fluctuation colour of the pop - 
ulation jBuzzomHl993l iBlakeslee. Vazdekis. fc Aiharll200lTl . 
or the mean of th e integrated and the fluctuation colours 
jTonrv et aljll997|) ]. With all these caveats, we have trans- 
for med the relat i on be tween Mfwow and (V — I c ), derived 
bv I Jensen et al] i2003fl from the 47 galaxies in their sample 
that show no signs of dust in the NIC2 field of view, into one 
between Mh and (V — I g ). To this end, we have relied on 
cq.|5| and on the transformation between Mh and Mfwow 
obtained in eq. 2 of Paper I. 11 For (J — K s ), we have sub- 
stituted 0.92 ± 0.06, which is the average integrated colour 
of the mentioned 47 galaxies, as obtained from the 2MASS 
Extended Source Catalog (XSC), via the GATOR catalog 
query page. This is the result: 

M„ = (-5.12 ± 0.09) + (4.9 ± 0.5) [(V - I g )o - 1.21]. (5) 

For the MC clusters, using the values for M h derived 
in Paper I and in lGonzalez. Liu, fc Bruzuall ll2005h . with the 
fluctuation magnitude for the supercluster type I once again 
corrected for extra extinction as in iJU the relation reads: 

M H = (-4.0 ± 0.5) + (4.8 ± 0.7) [(V - I g )o - 1.21]. (6) 

Even if tentative, since the near-IR observations of 
galaxies and star clusters have not been performed through 
the same filter, the conclusion is the same as the one drawn 
from the K s fluctuations: the slope for early-type galaxies 
and MC superclusters is the same within the errors; the 
MC clusters have fainter fluctuation magnitudes at a given 
(V — Ig) colour, as a result of their lower metallicity. It would 
certainly be worthwhile to obtain the data needed to per- 
form a fairer comparison of galaxies and MC star clusters in 
the H observing window in the future. Table 4 lists the coef- 
ficients of our fits to the MC supercluster data; the reduced 
chi-square (x 2 ) and the rms of the points (in magnitudes) 
after the fits are also included. 



5 SUMMARY & CONCLUSIONS 

In this paper we have presented the relation between abso- 
lute ifs-band SBF magnitude and (V — /) integrated colour 
for more than 180 MC star clusters. The newly reported re- 
sults extend to fluctuation magnitudes Mk b ~ -9 and optical 
colour (V — I) ~ -0.4 the linear relation observed already 
for early- type galaxies and spiral galaxy bulges in a more 
limited range of both fluctuation magnitudes and colour [-5 
> Mk. > -7, 1.05 <(V-I)< 1.25]. 

This empirical relation has been compared to the pre- 
dicted evolution of single age, single metallicity stellar pop- 
ulation properties in the Mr, vs . (V — I) diagram, based 
on bot h lBruzual fc Chariot] (|2003) and Mou hcine &: Lanconl 
(2003) isochrones. The predicted evolution and the observed 
sequence agree quite well over the ranges of Mk b magnitude 
and (V — 7) colour covered by the data [-5 > Mk b > -9, -0.4 

11 M h = M Fimw + (0.08 ± 0.07) - (0.24 ± 0.05) (J - K s + 
0.48). 
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Table 4. Fluctuation absolute magnitude vs. colour 



Ma b N rms \ 2 

M Kb -5.1 ± 0.4 4.0 ± 0.6 6 0.85 1.02 
M H -4.0 ± 0.5 4.8 ± 0.7 6 1.04 1.28 



Fits of the form M = a + b [ (V — I g )o — reference colour]; the number of objects used for the fit is 
tabulated as N. The resulting rms of the points (in magnitudes) after the fit and the reduced chi-square 
(x 2 ) are also listed. 



< (V— I) < 1.25]. With the exception of the extremely young 
supercluster Pre-SWB, and possibly the old and metal poor 
supercluster class VII, the remainin g superclus t ers lie on the 
linear correlation already found bv iLiu et alJ (|2002) . Such 
correlation is observed in the ranges -5 > Mk b > -9, 0.3 

< (V — I) < 1.25. The existence of a linear correlation im- 
plies that star clusters, early-type galaxies and spiral bulges 
represent an age sequence, where younger stellar popula- 
tions display brighter A's-band SBF magnitudes and bluer 
(V ~ I) colour. At the same time, the discernable offset be- 
tween galaxies and MC star clusters confirms that the Mk b 
vs. (V — /) plane may contribute to distinguish the effects 
of age and metallicity in intermediate-age and old stellar 
systems. One other suggestive result that emerges from the 
comparison between the galaxy sample and the models is 
that longer lifetimes of TP-AGB stars might partly explain 
galaxies with near-IR SBFs that are brighter than antici- 
pated, given their (V — /) colour. A preliminary compar- 
ison between the H 2MASS data of the MC star clusters 
and the sa mple of 47 ear l y-type galaxies and spiral bulges 
observed bv lJensen et al.H2003l) through the F1QQW filter 
leads to the same basic conclusions: galaxies and star clus- 
ters lie along correlations with the same slope, and there 
is a slight offset between the star cluster sample and the 
galaxies, caused by their different metallicities. 

We have found that results from star clusters in the 
MC (i.e., single stellar populations in Local Group irregular 
and consequently relatively metal-poor galaxies) agree with 
those of spiral bulges and early-type galaxies (i.e., compos- 
ite systems with higher metallicities than the MC), located 
not only in the Local Group, but also in mildly dense clus- 
ters of galaxies like Fornax and Virgo. The implication is 
that the relationship between M k s and (V — /) might be 
a fairly robust tool, rather insensitive to environment, at 
least in the local universe, for the study of ages and metal- 
licities of unresolved stellar populations; could provide ad- 
ditional constraints on star formation histories; and aid in 
the calibration of the i^-band SBFs for cosmological dis- 
tance measurements. In this regard -the determination of 
cosmological distances-, the sensitivity of near-IR SBFs to 
stellar ages and metallicities, and perhaps to the details of 
AGB evolution in intermediate-age populations, is a poten- 
tial caveat to bear in mind, albeit not necessarily in the 
local univers e or in view of today 's observing capabilities. 
For example, iFerreras et alJ lll999T) have analised early- type 
galaxies in Coma and 17 clusters at 0.3 < z < 0.9; their 
results imply that only galaxies smaller than 0.5 L„, at red- 
shifts z > 0.5, can be expected to harbour populations with 
metallicities below solar. 
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It is beyond the scope of the present paper to investi- 
gate in depth the role of stochastic effects in the integrated 
properties of stellar populations. However, for illustration 
purposes, we show in Fig. [I] the colour-magnitude diagrams 
(CMDs) of superclusters type IV (left panel) and VI (right 
panel). Fig.|K|is a copy of the right panel of Fig. [5] this time 
indicating the locations of individual clusters types IV (open 
triangles) and VI (open circles) in the Mk s vs. (V— I) plane. 
There are several facts worth noticing. Firstly, the colour 
magnitude diagrams, narrow at the top and wide at the 
bottom, resemble more those of SSPs with "normal" photo- 
metric errors than CMDs of composite populations (Bressan 
2005, private communication). This notwithstanding, and 
even when the relative offset between the general loci of the 
two SWB classes is clearly discernible, single clusters of each 
SWB type show a large scatter, 12 especially in fluctuation 
magnitude, which is determined by many fewer stars than 
the integrated colour. Taken together, the relatively nar- 
row RGB and AGB, and the very large scatter of individual 
clusters in the Mk s vs. (V — I) plane support the interpre- 
tation that stochastic effects might be more important than 
actual population variations in determining the integrated 
properties of single clusters within each SWB class. If this is 
the case, the construction of superclusters is an appropriate 
strategy to simulate more massive SSPs. 

This paper has been typeset from a TjrjX/ DTfrjX file prepared 
by the author. 



On the basis of this plot alone, it could be argued that per- 
haps SL 663 , at (V — /) ~ 1.7, does not belong in class SWB 
IV. Removing it leaves the fluctuation magnitude and integrated 
colour of supercluster type IV the same, within the quoted errors. 
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Figure 4. Colour-magnitude diagrams of MC superclusters types SWB IV (left panel) and SWB VI right panel). Stars within 60 arcsec 
from the center (at the distance of the LMC). Average photometric errors are 0.04 mag in brightness and 0.02 mag in color for sources 
with K s ^ 13; 0.06 and 0.03 mag for stars with 13 < K B ^ 14; and 0.13 and 0.07 mag (about the size of the dots) for sources with 14 
< K s < 15. 




Figure 5. Mk, versus (V — I g ) colour of individual clusters classes SWB IV (open triangles) and SWB VI (open circles); all other 
symbols as in right panel of Fig. 



